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A two-phase model for the Narrow Line Region of NGC4151 
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Abstract. NGC4151 is one of the brightest and best-studied Seyfert galaxies. Here, we present a two-phase model 
of the narrow-line region (NLR) of NGC4151. This study is motivated by (i) the fact that the X-ray spectrum of 
NGC4151 is among the flattest known for Seyferts, and (ii) the recent Chandra detection of an X-ray narrow-line 
region in this galaxy. X-ray spectra as flat as that of NGC4151 (r x ~ —1.5) are expected to favor the presence 
of two gas phases in pressure equilibrium (Krolik, McKee & Tarter 1981). In the present study, we show that 
a pronounced two-phase equilibrium develops in the extended emission-line region if we use the observed multi- 
wavelength spectrum of NGC4151 to ionize the clouds. The material is stable to isobaric perturbations over a 
wide range of temperatures. We therefore propose that such a condition has arisen in the NLR of NGC 4151, and 
that it explains the detection of hot, extended X-ray gas which we identify as the NLR-cloud confining medium. 
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1. Introduction 

1.1. NGC4151 

NGC 4151 is the optically brightest Seyfert galaxy of 
type 1.5. It is very nearby (r ~ 20 Mpc for Hq=50 
km/s/Mpc) and is therefore excellently suited to study in 
detail the physical processes in the nuclei of active galaxies 
( AGN) . The extended narrow-emission-line region (NLR) 
of NGC 4151 has a cone-like geometry (Schulz 1988, 1990, 
Evans et al. 1993, Boksenberg et al. 1995), interpreted as 
'radiation cone' (Schulz 1988) caused by anisotropic illu- 
mination of the NLR. In the context of the unified model 
of AGN (Antonucci 1993) such a radiation geometry arises 
from the partial obscuration of a molecular torus or the 
anisotropic emission of the accretion disk. 

NGC 4151 is the most intensely studied AGN across 
the electromagnetic spectrum (see Schulz 1995 and Ulrich 
2000 for detailed reviews). The X-ray spectrum of the 
galaxy is complex in the soft band (e.g., Weaver et al. 
1994) and extends to very high energies (e.g., Perotti et 
al. 1981, Maisack et al. 1993). The powerlaw spectrum 
has a mean photon index T x ~ —1.5 and is significantly 
flatter than that of other Seyfert galaxies (<T X >~ —1.9). 
Using the Chandra HETG spectrometer, Ogle et al. (2000) 
recently reported the detection of a high-temperature, 
narrow-line, X-ray emitting plasma in NGC 4151, confirm- 
ing earlier Einstein (Elvis et al. 1983) and ROSAT (Morse 
et al. 1995) observations of extended X-ray emission in 
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this galaxy. The X-ray gas is spatially coincident with the 
NLR and extended narrow- line region (ENLR). A hot in- 
tercloud medium was one of several possible explanations 
mentioned repeatedly (e.g., Elvis et al. 1983, 1990, Morse 
et al. 1995, Ogle et al. 2000). It is therefore important 
to examine if, and under which conditions, a two-phase 
equilibrium of the NLR and ENLR clouds is possible. 



1.2. Multi-phase cloud models, and motivation of the 
present study 

Different mechanisms for the confinement of the emission- 
line clouds in AGN have been studied over the last few 
decades. The model originally investigated in detail by 
Krolik, McKee & Tarter (1981; KMT hereafter) is a 
two-phase model, consisting of cold line-emitting clouds 
(T ~ 10 4 K) in pressure balance with a hot inter-cloud 
medium (T ~ 10 8 K). KMT assumed a relatively flat X- 
ray spectrum of their input continuum illuminating the 
clouds. Later it turned out that Seyfert galaxies typically 
show steeper X-ray spectra and a pressure balance be- 
tween a cold, photoionization heated, and a hot, Compton 
heated phase then no longer exists (e.g., Fabian et al. 
1986). More recent studies have demonstrated the pres- 
ence of an additional stable region of intermediate tem- 
perature (Reynolds & Fabian 1995, Komossa & Fink 1997, 
Komossa 2001), and the recent discovery of warm ab- 
sorbers located in that intermediate region has revived 
the interest in multi-phase cloud models. 
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The aim of this paper, however, is not a general dis- 
cussion of the viability of variants of the KMT model to 
explain the BLR/NLR-confincmcnt in Seyfert galaxies in 
general. Instead, we address the following question: will 
a KMT-like two-phase equilibrium develop, if conditions 
are favorable in individual galaxies ? 

The observations of extended X-ray emission in 
NGC4151 spatially coincident with the NLR, together 
with the fact that the relatively hard X-ray spectrum of 
NGC4151 is favorable for the development of a two-phase 
equilibrium, motivated the present study. 
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Fig. 1. Average radio to gamma-ray continuum spectrum 
of NGC4151, composed from more than 300 individual 
measurements collected from the literature. The unob- 
served EUV part of the spectrum is shown as dotted line. 
The hard X-ray spectrum was conservatively assumed to 
show a break, whereas most measurements indicated an 
extent of the spectrum to much higher energies (see Sect. 
3.1.1). 



2. Model calculations and results 

The stability of photoionized gas clouds to isobaric pertur- 
bations can be examined by studying the behavior of cloud 
temperature T as a function of pressure. If the tempera- 
ture is multi-valued for constant pressure, and the gra- 
dient of the equilibrium curve is positive, several phases 
may exist in pressure balance. 

The photoionization calculations were carried out with 
the code Cloudy (Ferland 1993) under the following as- 
sumptions: We adopted solar metal abundances according 
to Grevesse & Anders (1989). The gas density was as- 
sumed to be constant (logriH = 3), as was the column 
density (logiVn = 21.5). Test calculations with differ- 
ent densities and abundances were performed but they 
do not alter the results discussed here. The clouds were 
illuminated by the continuum of a central point-like en- 
ergy source. The spectral energy distribution (SED) used 
for modeling corresponds to the average multi-wavelength 



continuum of NGC4151 (Fig. 1), collected from the lit- 
erature and determined from more than 300 individual 
data points (Komossa 1994, Schulz & Komossa 1993, and 
references therein) . The adopted X-ray spectrum is partic- 
ularly important in deriving the shape of the equilibrium 
curve (T versus U/T) for high temperatures. We used the 
mean observed powerlaw index, T x = —1.5 (e.g., Perotti et 
al. 1990), extending up to 100 keV, where we assumed the 
spectrum to break steeply into the gamma-ray regime (see 
Sect. 3.1.1 for further details and references). This mean 
SED is appropriate for modeling, because the recombina- 
tion time scale for gas of NLR densities is long, and the 
NLR thus sees an average continuum. 

The cloud temperature as a function of the ionization 
parameter was extracted after each model run. The ion- 
ization parameter is defined as 

oo 

U = — % , with Q = And 2 f p-du , (1) 

47rr^nfjc J hv 

where rin is the gas density, r the distance of the gas cloud 
from the continuum source, Q the number rate of photons 
above the Lyman limit vq, and d the distance between 
observer and galaxy. 

Results are shown in Fig. 2. We find that a clear two- 
phase equilibrium between a low-temperature and a high- 
temperature component develops over a range of ioniza- 
tion parameters. We have verified that our results are 
independent of metal abundances^], gas density, and the 
amount of absorbed radiofrequency radiation. 

In a next step, we added a black-body component to 
the UV part of the spectrum. A SED with a huge EUV 
excess which dominates Q (Qhb — 0.8Qtoti) was invoked 
by Penston et al. (1990) using the ionization-parameter- 
sensitive emission-line ratio [OII]/[OIII]. This SED would 
significantly narrow down the region where a two-phase 
equilibrium is possible (Fig. 2). Since Schulz & Komossa 
(1993) showed that the [OII]/[OIII]-ratio overestimates 
the ionization parameter (thus Q) if the emission-line re- 
gion is inhomogeneous in density, and since recent studies 
did not favor a giant EUV excess in this galaxy (Alexander 
et al. 1999), we do not discuss this SED further. 

For comparison, we have plotted in Fig. 2 the phase 
diagram obtained using the average SED of a sample 
of Seyfert galaxies (Komossa & Schulz 1997). This SED 
is characterized by a much steeper X-ray spectrum (L x 
=—1.9) and shows that the temperature of the upper 
branch has reduced and the overlap between upper an 
lower branch has disappeared. 

3. Discussion 



For supersolar metal abundances the small region at T ~ 
10 6 K which allows a three-phase equilibrium becomes more 
prominent. However, photoionization modeling of the ENLR 
emission-line ratios of NGC4151 (Schulz & Komossa 1993) fa- 
vored depleted metallicity. 
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Fig. 2. Equilibrium gas temperature T as a function of 
U/T, obtained using the observed SED of NGC 4151 (thick 
solid line). The open circle marks a value of log U — —2.0. 
Two curves are shown for comparison. Dashed line: SED 
with an additional strong EUV-bump component (Qbb = 
0.8 Qtoti! see text); dotted line: mean Seyfert SED with T x 
=-1.9. 



3.1. Model assumptions 

Before we discuss consequences of the models considered 
here, we first give some cautious comments on the gen- 
eral model assumptions. The assumptions of photoioniza- 
tion equilibrium and thermal balance have to be kept in 
mind, as already discussed by KMT. For instance, the gas 
cooling time scale should be less than the dynamical time 
scale (Morse et al. 1995 estimate that the time scales are 
comparable in the ENLR of NGC 4151). The assumption 
of photoionization equilibrium is problematic if intrinsic 
source variability and recombination time scale are of the 
same order. Given the low mean density of the NLR and 
ENLR of NGC 4151 (logra ~ 2-3) the recombination time 
scale, i r oc,H = a _1 (T e )n~ 1 , is on the order of 1.3 10 2-3 
yrs, and correspondingly longer for even lower density. 
This is long compared to the rapid X-ray variability of 
this source. An average continuum is thus appropriate for 
modeling. Heating agents in addition to the central con- 
tinuum source (like cloud friction) were neglected. KMT 
showed that their influence is generally weak. However, 
radiation pressure by trapped line photons could signifi- 
cantly increase cloud pressure (Elitzur & Ferland 1986). 

Apart from these general considerations, the influ- 
ence of the shape of the EUV-X-ray continuum specific 
to NGC 4151 on the results is particularly important. 

3.1.1. X-ray continuum shape 

Relevant for the present study is the shape of the X-ray 
spectrum (its steepness, and the position of the break). 
We used an input spectrum with T x = —1.5, the mean ob- 



served spectral index (e.g., Perotti et al. 1990)0, extending 
up to 100 keV. At 100 keV the spectrum shows a break 
(Maisack et al. 1993) and continues with a significantly 
steeper powerlaw. It is interesting to note that early X-ray 
observations of NGC 4151 detected the source out to much 
higher energies (e.g., Schonfelder 1978, Perotti et al. 1981, 
Perotti et al. 1991), whereas some recent observations fa- 
vored a break at 50-100 keV (Apal'kov et al. 1992, Maisack 
et al. 1993). Were the early detections of the source based 
on measurement errors ? We do not believe so, because a 
few early observations did show a break in the spectrum 
of NGC 4151 (Baity et al. 1984, Bassani et al. 1986). In 
particular, Baity et al. presented several measurements, 
of which only one showed a spectral break. Since taken 
with the same instrument, but at different epochs, mea- 
surement errors are unlikely, and it is more plausible that 
NGC 4151 underwent real spectral variability. 

To be on the 'safe side', we conservatively assumed the 
spectrum to break, at inbreak = 100 keV. Had we adopted 
a spectrum extending up into the Gamma-regime (MeV 
energies) the two-phase equilibrium would have become 
even more pronounced, due to the increased importance 
of Compton heating. 

3.1.2. EUV spectrum 

Secondly, the EUV spectral shape is of importance for the 
shape of the equilibrium curve. Using a continuum with a 
strong EUV excess reduces the Compton temperature of 
the upper branch and thus narrows down the range were a 
two-phase equilibrium is possible (Fig. 2; see Fabian et al. 
1986, Komossa & Fink 1997, Komossa & Meerschweinchen 
2000 for details). 

Many different EUV continuum shapes of NGC 4151, 
and the ENLR emission-line response to them, were exam- 
ined by Schulz & Komossa (1993). They concluded that 
the observed ENLR line-ratios can be explained by two 
very different models: (i) a SED with a huge EUV bump 
as suggested earlier by Penston et al. (1990) based on Q es- 
timates, and (ii) a powerlaw-like EUV continuum of much 
lower Q, close to the lower limit set by the Hf3 constraint 
(and a density-mixed, partly matter-bounded ENLR, to 
produce the same observed emission-line ratios). A con- 
tinuum of the latter form (ii) was favored by Alexander et 
al. (1999), and adopted by Kraemer et al. (2000), and we 
also use it in the present analysis. 



2 for comparison: the following photon indices (either means 
over a collection of observations, or single measurements) of the 
hard X-ray spectrum of NGC 4151 were reported in the litera- 
ture: Beall et al. (1981): T x = -1.4, Baity et al. (1984): -1.6, 
Warwick et al. (1989): -1.54, Maisack and Yaqoob (1991): 
-1.37, Ogle et al. (2000): -1.4; for a collection of T x mea- 
surements of NGC 4151 prior to 1989 see Tab. 2 of Baity et al. 
(1984), and Tab. 2 of Warwick et al. (1989) 
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3.2. Two-phase NLR of NGC 4151 

Over the years, a number of different models for the 
emission-line regions of AGN have been suggested. Most 
concentrated on the BLR. The broad-line clouds might 
be confined by mechanisms other than thermal pressure 
(e.g. magnetic fields could play a role; Rees 1987), or they 
might not be confined at all. Confinement is unnecessary if 
the BLR clouds are continuously produced (e.g., Perry & 
Dyson 1985, Murray & Chiang 1995). Generally, all sug- 
gested models including variants of the KMT approach, 
still suffer from some unsolved problems. 

In the present study, it is not our aim to find the ulti- 
mate BLR/NLR model relevant for all types of AGN un- 
der any conditions. We rather address the question: will 
a KMT-like two-phase equilibrium develop in individual 
galaxies if conditions are favorable ? The previous obser- 
vations of NGC 4151 are highly suggestive that this is in- 
deed the case. Firstly, the X-ray spectrum of NGC 4151 
is very hard. Secondly, extended X-ray emission, which 
could play the role of the NLR/ENLR confining medium, 
has already been detected. The Chandra observation re- 
vealed similar kinematics of optical and X-ray NLR which 
suggests that they are closely linked (Ogle et al. 2000). 

A number of explanations for the hot extended X-ray 
emission of NGC 4151 have been discussed (Elvis et al. 
1983, 1990, Heckman & Balick 1983, Morse et al. 1995): 

(i) emission from shocks due to NLR-cloud - ISM inter- 
action. However, this would imply a very small covering 
fraction of the X-ray gas, and would not work for the 
ENLR, since cloud velocities (Schulz 1988) are too small. 

(ii) Scattered nuclear emission was considered a possibil- 
ity, but Morse et al. (1995) presented estimates which 
make this scenario unlikely, (iii) Thermal emission from 
a hot confining medium has been previously mentioned; 
and we have shown here that a two-phase medium devel- 
ops if illuminated by the observed SED of NGC 4151. 

Some important consistency checks: The mean den- 
sity of the ENLR, as inferred from the density-sensitive 
sulphur emission-line ratio, is logn = 2.3 (Schulz 1988, 
Penston et al. 1990). There are indications that the density 
increases inwards (e.g., Schulz 1995). For the estimates be- 
low, we adopt the generally accepted density scaling law of 
the NLR, office \, with logn co id = 2.3 at 1.5 kpc (corre- 
sponding to lognhot = —1.9; see Fig. 2) and logn co id = 6.3 
at 15 kpc (logrihot = 2.1). What is the column density of 
the hot inter-cloud medium ? The ENLR of NGC 4151 ex- 
tends out to ~20"(=1.9 kpc), similar to the distance to 
which the X-ray emission has been traced (1.5 kpc, Morse 
et al. 1995) in the south-west cone. We then expect a col- 
umn density of log ATh — 6 10 21 cm -2 of the hot gas. The 
medium is therefore optically thin to Thomson scattering 
and would not produce any observable effects which have 
not been detected. The observed X-ray luminosity of the 
extended component amounts to ~10 41 erg/s (Elvis et al. 
1983 : 5 1 40 erg/s, Morse et al. 1995: 3 10 41 erg/s). This 
compares to an X-ray luminosity of 1 10 41 erg/s predicted 
by our model, assuming an NLR/ENLR extent and den- 



sity law as given above, and a total filling factor of the 
gas that corresponds to a cone opening angle of 75° (e.g., 
Evans et al. 1993). We note that this comparison should 
be regarded an order-of-magnitude estimate, due to un- 
certainties in both, observations and modelling assump- 
tions. The hot gas phase is highly ionized and we predict 
the dominance of lines from the ions of, e.g., CVI, NVII, 
OVIII, SiXIII-XIV, and FeXXIV-XXVI. Ionization stages 
up to FeXXV have indeed been detected by Chandra. 

4. Summary and conclusions 

The close proximity and brightness of NGC 4151 make it a 
uniquely valuable object to explore the interaction of the 
active nucleus with its immediate environment. 

We have assessed the influence of the observed flat X- 
ray continuum of NGC 4151 on the thermal stability of 
NLR/ENLR clouds and the possibility of cloud confine- 
ment by a hot inter-cloud medium. We find that a two- 
phase medium develops and we suggest to identify the 
hot phase with the observed extended X-ray emission in 
NGC 4151, which is co-spatial with the NLR/ENLR. 

We estimate a low density of the extended gas which 
is shown to be able to confine the total NLR and ENLR 
without becoming Thomson-thick. 

The results presented here are relevant for other 
nearby Seyfert galaxies with relatively hard X-ray spectra 
(e.g., NGC 3227: T x k —1.5). Such galaxies are therefore 
important targets for future Chandra and XMM-Newton 
observations. 
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